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ABSTRACT 

We present an analysis of a 72 ks Chandra observation of the double cluster Abell 1644 (z=0.047). 
The X-ray temperatures indicate the masses are M 500 = 2.6 ± 0.4 x 10 14 /i _:L M Q for the northern 
subcluster and M500 = 3.1 ± 0.4 x 10 14 Ii~ 1 Mq for the southern, main cluster. We identify a sharp 
edge in the radial X-ray surface brightness of the main cluster, which we find to be a cold front, with a 
jump in temperature of a factor of ~3. This edge possesses a spiral morphology characteristic of core 
gas sloshing around the cluster potential minimum. We present observational evidence, supported 
by hydrodynamic simulations, that the northern subcluster is the object which initiated the core gas 
sloshing in the main cluster at least 700 Myr ago. We discuss reheating of the main cluster's core 
gas via two mechanisms brought about by the sloshing gas: first, the release of gravitational potential 
energy gained by the core's displacement from the potential minimum, and second, a dredging inwards 
of the outer, higher entropy cluster gas along finger-shaped streams. We find the available gravitational 
potential energy is small compared to the energy released by the cooling gas in the core. 
Subject headings: galaxies: clusters: general; galaxies: clusters: individual (Abell 1644); galaxies: 
interactions; Xrays: galaxies: clusters; galaxies: interactions 



1. INTRODUCTION 

Double and multiple cluster systems arise as a natu- 
ral consequence of hierarchical structure formation. Be- 
ing acutely sensitive to the gas density in clusters, X- 
ray observations provide a powerful method for unam- 
biguously identifying cluster systems and substructure 
within them. With early X-ray telescopes such as Ein- 
stein, ROSAT and ASCA, 60% of clusters were catego- 
rized as "relaxed" , based on the smooth, mostly circular 
flux isophotes centered on their cores (e.g. Jones & For- 
man 1999). It is only with the most recent generation of 
X-ray telescopes (XMM-Newton and, to a greater degree, 
Chandra) that the angular resolution has increased suf- 
ficiently to examine small scale features near the cluster 
cores. Many clusters previously considered relaxed, upon 
closer inspection, reveal highly disturbed cores. The dis- 
turbances include nearly ubiquitous AGN-blown bubbles 
in the intracluster gas (c.f. Churazov et al. 2000, Fabian 
et al. 2000, Nulsen et al. 2002, and review in McNamara 
&: Nulsen 2007) and cold fronts (review by Markevitch 
& Vikhlinin 2007, hereafter MV07). Markevitch et al. 
(2003) found that more than 2/3 of the clusters classified 
as cool core systems in a flux limited sample by Peres et 
al. (1998) exhibit edges in the X-ray surface brightness 
profiles around their cores when viewed with Chandra. 
X-ray spectra show that most of these edges are contact 
discontinuities, or "cold fronts" (MV07), though occa- 
sionally they are caused by shock fronts (e.g. Hydra A, 
in Nulsen et al. 2005). 

Cold fronts indicate gas motion in a cluster and may 
be located at the boundary of an infalling subcluster, 
as in the merger in A3667 (Vikhlinin et al. 2001) or in 
1E0657-56 (Markevitch et al. 2002). They may also arise 
in the core of the primary cluster during the passage of 
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a subcluster, when the core gas starts "sloshing" in re- 
sponse to the gravitational disturbance (Markevitch et 
al. 2001; Churazov et al. 2003; Ascasibar & Markevitch 
2006, hereafter AM06). We find direct evidence in sup- 
port of this gas sloshing in the core of Abell 1644, which 
we explore in this paper. 

2. CORE GAS OSCILLATION (SLOSHING) 

"Sloshing" of the dense cluster core gas was first pro- 
posed by Markevitch et al. (2001) to explain a cold front 
observed in the core of a relaxed cluster, Abell 1795. 
This process was studied via hydrodynamic simulations 
in AM06 and reviewed in detail in MV07. A similar pro- 
cess to explain these features, though slightly different in 
interpretation (MV07), was proposed by Tittley & Hen- 
ricksen (2005). 

From the above works, it was found that core sloshing 
in otherwise undisturbed clusters arises when the cluster 
undergoes a perturbation to its gravitational potential 
caused by another infalling group or cluster, even if that 
infalling cluster is gasless and the resulting large scale 
disturbance is small. The sloshing occurs as a result of 
the gas core lagging behind the cluster the potential min- 
imum, as they both move towards the perturbing object. 
As the gas core falls back onto the potential minimum, 
it overshoots it and begins to oscillate. With each os- 
cillation, the gas core is moving against its own trailing 
gas, producing an "edge" in the X-ray brightness which 
expands out from the cluster. This sequence of events is 
described in more detail in AM06. 

The continued oscillation of the core gas about the po- 
tential minimum produces a succession of radially prop- 
agating cold fronts, manifested as concentric edges in the 
surface brightness distribution. These fronts may form 
a spiral structure when the sloshing direction is near the 
plane of the sky and the merger has a non-zero angu- 
lar momentum (AM06). A recent study by Lagana et al. 
(2009) examining substructure in a sample of X-ray clus- 
ters has found such a spiral pattern around the cores of 
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those clusters. Here, we see exactly this spiral structure 
in the core of the main cluster in Abell 1644. 

In § 3, we describe the reduction of the Chandra data 
on A1644 and extraction of X-ray spectra. In § 4 we list 
the various gas quantities from the fitted models. In § 5 
we discuss implications of these measurements, and our 
conclusion is in § 6. 

3. PREVIOUS AND NEW OBSERVATIONS OF A1644 

We present a 72 ks Chandra ACIS-I observation (Ob- 
sIDs 7922 and 2204) of the cluster A1644 (z=0.047), a 
rich binary cluster located ~3 Mpc from the Shapley Su- 
per cluster (De Filippis et al. 2005), likely lying along 
one of the Shapley filaments (see Figure 1 and Figure 2). 

The most complete sample of cluster members in 
A1644 contains 141 galaxies (Tustin et al. 2001). They 
find a lack of evidence for significant substructure in the 
galaxy distribution and so do not differentiate between 
the main and sub-cluster but instead give a mean redshift 
cz = 14,295 ± 93 km s _1 and <7„=1108 km s^ 1 for the 
entire cluster. Although there is disagreement in the two 
previous redshift studies (Dressier & Schectman 1988; 
Tustin et al. 2001) as to the presence of substructure, 
A1644 has a clear bimodal distribution in X-rays, first 
observed with the Einstein X-ray Observatory (Jones & 
Forman 1984). 

Subsequent XMM-Newton observations (Reiprich et al. 
2004) revealed that the two X-ray clusters (a main- and 
sub-cluster) are connected by a warm (4-6 keV) , tenuous 
intracluster medium. The two X-ray peaks coincide with 
two of the brightest galaxies in the region. The south- 
ern, brighter X-ray peak lies on the brightest galaxy in 
the cluster, a cD with M#= -25.84, while the northern, 
fainter X-ray peak lies on the 5th brightest galaxy, an el- 
liptical, with M H = -24.74 (Tustin et al. 2001). Thus we 
assume that the X-ray peaks and the galaxies all lie at lo- 
cal minima in the cluster gravitational potential. Based 
on the H-band magnitudes (Tustin et al. 2001), and the 
X-ray data, we henceforward refer to the SW cluster as 
the "main" cluster or A1644-S and the NE cluster as the 
"subcluster" or A1644-N. 

3.1. ACIS-I Data Reduction 

A1644 was observed with ACIS-I, using chips 10-3, S2, 
S3 in two observations (ObsIDs 7922 and 2204) with to- 
tal exposure time of 72 ks. We reprocessed the level 1 
event files as described in Vikhlinin et al. (2005). The 
background flare rejection was done using the ratio of 
the 2.5-7 keV and 9.0-12 keV bands (Hickox & Markc- 
vitch 2006) over the CCDs 10-13 and S2, excluding bright 
source and cluster emission. This resulted in a total ^70 
ks clean exposure, on which we performed further spec- 
tral and image processing tasks described below. Among 
the products created in this initial reduction are back- 
ground images and exposure maps for several interesting 
energy bands, which incorporate the telescope aspect so- 
lution, vignetting, energy dependent effective area, and 
position and energy dependent detector efficiency. 

Proper estimation of the X-ray background is critical 
for spectral analysis of the faint extended emission from 
clusters. We used the blank sky datasets from the Chan- 
dra calibration database (CALDB). The aspect solution 
for our observation was applied to the blank sky expo- 
sures, remapping the background events to our source 



coordinate frame. We then performed an overall back- 
ground normalization by taking the ratio of the count 
rates at high energies (9-12 keV, where the telescope's 
effective area drops dramatically) in our source images 
to that in our background images. For each of our spec- 
tral extraction regions, an identical region on the detec- 
tor was chosen from the blank sky image to provide a 
background estimate. Associated ARFs and RMFs were 
created using CI AO v4.1. 

3.2. X-ray Imaging and Photometry 

From the Chandra event file, we created an image in 
the 0.5-2.5 keV band in order to maximize detector sen- 
sitivity and contrast of the soft cluster emission against 
the detector background. We then subtracted the back- 
ground image and divided the result by the correspond- 
ing exposure map. This exposure-corrected image was 
used for all subsequent position and imaging analyses 
and is shown in Figure 1 and Figure 2 (left). We de- 
tected point sources using the wavelet detection algo- 
rithm wvdecomp outlined in Vikhlinin et al. (1998). The 
detected point source regions were checked by eye and 
then excluded from all subsequent analyses. 

3.3. Spiral Structure in the Core of A 1644 S 

Chandra 's A CIS-I field of view is large enough to im- 
age both A1644-N and A1644-S on four CCDs (10-13), 
allowing a simultaneous view of both the ICM emission 
connecting the two subclusters and a high resolution pic- 
ture of their cores. In A1644-S, we see a spiral pattern in 
the core X-ray emission, defined by an edge just to the 
north of the brightness peak and extending over more 
than 180° in azimuth (see Figure 3 left). To examine the 
physical properties of the gas on either side of this promi- 
nent surface brightness edge near the core of A1644-S, we 
extract spectra in regions which trace this edge. In an 
effort to trace the edge with circular regions, which we 
require for the application of local spherical symmentry 
in §4.3, and to maximize our signal in each region, we 
choose the two sectors shown in Figure 4 (left). 

3.4. Spectral Extraction and Fitting 

Spectral extraction was performed using the CIAO 3 
tool specextract. We defined two sectors over which to 
select our regions, as the projected surface brightness 
edge could be adequately traced by a circular region over 
~60 degrees in azimuth. Within each sector, spectra 
were extracted in circular segments at varying radii from 
the cluster center. We also examine a surface brightness 
edge in the smaller subcluster to the north, over a sector 
centered on the X-ray peak where the isophotes appear 
compressed (see Figure 4, right). 

For fitting, spectra in each region was binned with a 
minimum of 30 cts per bin in the inner regions (SQ2, 
SQ3 and NQ1), and 50 cts per bin in all other regions. 
For each spectrum, we fit a thermal MEKAL model con- 
volved with an intergalactic absorption model (xsmekal * 
xszwabs). The cluster's redshift was fixed at the median 
redshift (z=0.047) obtained from the spectroscopic sam- 
ple in Tustin et al. (2001) with an absorbing hydrogen 

3 All spectral extraction and spectral analysis were performed 
using CIAO v4.1, the CALDB v4, and Sherpa. 
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Fig. 1. — Combined 70ks Chandra ACIS-I image (obsid 2206 and 7922) over the energy range 0.5-2.5 keV. We see the two bright cores of 
A1644-N and A1644-S. Both the main and subcluster's emission is preferentially extended towards the other, indicating previous interaction 
between the two systems. We also see the spiral pattern of the gas sloshing in A1644-S. The image has been background subtracted, exposure 
corrected, binned by a factor of 2 and smoothed in ds9 using a Gaussian kernel with radius 4 pixel (~ 4"). 
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column density N# = 5.3 x 10 20 cm -2 (Dickey & Lock- 
man 1990) at the cluster position. We account for the 
gas lying in projection along the line of sight by using a 
variation of the algorithms projct for XSpec and depro- 
ject for Sherpa. The basic algorithm is to take a series of 
radially extracted spectra, fit a model to the outermost 
bin first and then move inward. For each inward step, an 
additional model is added to the fit, whose gas temper- 
ature (and optionally abundance) is fixed at the values 
obtained from the outer bins. The difference between 
these and our method lies in our separate derivation of 
the 3D gas density rather than leaving the density as a 
free parameter in the model fit. The density for each ra- 
dial bin is calculated as described in §4.3 and the ratio of 
the emission measure in the projected region to that in 
the emission region is used as a normalization constraint 
for each of the projected regions. 

The deprojection of this system is complicated by 
the extended cluster emission that is preferentially dis- 
tributed along the merger axis and largely fills the ACIS- 
I camera (see Figure 1 and Figure 2, right). So, simply 
extending our spectral extraction regions (particularly 
SQ2) out to large radii will inaccurately predict the emis- 
sion along the line of sight. Instead, we examine sectors 
orthogonal to the merger axis of the cluster with radii 
out to the detector edge, which are likely to be more 
representative of the line of sight gas distribution. We 
then fit a density model as described in §4.3 and used 
this to calculate the emission measure for the outermost 
regions in the deprojection (labelled Nrjeproj and Soeproj 
in Figure 2). The model temperatures are listed in Ta- 
ble 1, plotted in Figure 5 and Figure 6 for each region in 
Figure 4, and are discussed below. 

4. CLUSTER MEASUREMENTS 
4.1. Main and Subcluster Data 

In Table 2, we show the best fits for several large re- 
gions around the cluster. In the circular regions sur- 
rounding, but excluding, each of the cool cores, we find 
kT=4.62t°; 2 ^ keV for A1644-N and kT=5.10±0.14 keV 
for A1644-S. These regions are shown in Figure 2 (left). 
We compare the Chandra and XMM-Newton tempera- 
tures using the portion of the XMM-Newton region re- 
ported by Reiprich et al. (2004) that is covered in the 
smaller Chandra FOV. Here, we obtain kT=4.42 ± 0.05 
keV and abundance 0.44 ± 0.03, as compared to their 
kT=3.83 ± 0.06 keV and abundance 0.32+^;° 2 . Abun- 
dances are reported with respect to solar using the mod- 
els of Anders & Grevesse (1989). 

Using the Msoo-Tx relation from Vikhlinin et al. 
(2009), we find M 500 = 2.6±0.4x lO^h^M® for A1644- 
N and M 50 o = 3.1 ± 0.4 x IO^H^Mq for A1644-S. The 
latter mass is consistent with the mass derived for A1644- 
S from a caustic analysis in Tustin et al. (2001), where 
they find M 500 ~ 4.6 ± 0.7 x 1O 14 M 

4.2. Temperature Maps 

In Figure 7 we present temperature and temperature 
error maps for A1644-S and A1644-N. Temperature maps 
can reveal interesting structures and inhomogeneities in 
the cluster gas which might otherwise be overlooked. We 
use the method of Markevitch et al. (2000), in which a 
series of 6 images in non-overlapping energy bands (0.5- 



TABLE 1 
Radial Profiles 



Region 


(") 


(") 


K 1 

(keV) 


X ( d °t) 


^Deproj 


515.0 
350.0 
200.0 


700.0 
515.0 
350.0 


r ,-9+0.39 

a oc+0.39 
4-85-0.32 
5 D4+ ' 40 


0.97(269) 
0.75(263) 
0.76(229) 


SQ2 


120.0 
48.0 
12.1 
4.9 


200.0 
120.0 
48.0 
12.1 


5 91+O.37 
° -0.36 
5 36+ - 51 
c 7/ . + 1.23 

i.86tr 2 ? 


0.71(183) 
0.77(121) 
0.64(36) 
0.75(11) 


SQ3 


100.0 
32.6 
19.0 
9.0 


200.0 
100.0 
32.6 
19.0 


5.471^ 

c 77 +0.80 
°- ' '-0.57 
1 Q2+ - 13 

1 7c;+0-23 
L - ,o -0.14 


0.68(193) 
0.71(92) 
0.95(37) 
1.90(23) 


■N Deproj 


450.0 
250.0 
100.0 


600.0 
450.0 
250.0 


6.99tS;$ 
6.06lJ;i? 


0.75(199) 
0.83(241) 
0.62(120) 


NQl 


60.0 
19.7 
10.3 
4.9 


100.0 
60.0 
19.7 
10.3 


^• oo _0.68 
9 oc+0.67 
2-8t>_ 45 

1 40+ - 11 
1 -^ u -o.io 
95+0-08 


1.11(51) 
0.60(39) 
0.96(12) 
0.77(7) 



Note. — Table 1: MEKAL model temper- 
atures for the regions shown in Figure 2 (left) 
and Figure 4. ( col 1 ) Spectral extraction region 
identifier for regions shown in Figure 4 with the 
exception of the deprojection regions, which are 
shown in Figure 2 (left). (col 2 and 3) Inner 
and outer radii of sectors, (col 4) Deprojected 
MEKAL model temperature. Errors on kT are 
the 68% confidence interval, A\ =1, for one inter- 
esting parameter), (col 5) Reduced \ 2 f° r model 
and number of degrees of freedom. 



TABLE 2 

Spectral Fits for Large Regions of Interest 



Region 


kT 
(keV) 


Abundance 


X 2 (dof) 


Entire cluster 
A1644-N* 
A1644-S* 


4.42±0.05 
4 62+ ' 24 
5.10±0.14 


0.44±0.03 
0.44±0.08 
29+ 05 


1.66(804) 
1.66(403) 
1.06(576) 


SROI1 


3 46 +(, ' 2l) 


0.44 


1.10(63) 



Note. — Table 2: MEKAL model parameters for 
several large regions: the portion of the cluster region 
described in Reiprich et al. (2004), two circular re- 
gions surrounding each cluster's X-ray peak* (shown 
in Figure 2, left), and SROI1 (shown in Figure 4, left). 
(col 2) Temperature, (col 3) abundance, and (col 4) 
reduced x 2 from the MEKAL model fits. Errors are 
the 68% confidence interval on two interesting param- 
eters, Ax 2 =2.3, with other parameters free. Abun- 
dance for SROI1 was fixed at the cluster mean. *Ex- 
cluding the cool cores (r=56 kpc for N and r=97 kpc 
for S) 

1.5, 1.5-2.3, 2.3-3.5, 3.5-4.5, 4.5-6.5, 6.5-9.5 keV) are cre- 
ated along with corresponding exposure maps. Each sur- 
face brightness image is then smoothed identically with 
a variable length scale Gaussian kernel, such that the 
smoothing scale at each pixel decreases roughly as the 
square root of the number of counts in that pixel (out to 
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Fig. 2. — (left) Same as Figure 1 with A1644-N and A1644-S cluster regions, whose fits are in Table 2, shown as circles with the cores 
excluded along with outer deprojection spectral regions shown as sectors, (right) corresponding R-band optical image with X-ray surface 
brightness contours overlayed, showing the peaks of the X-ray emission lying on top of two bright elliptical galaxies. A low surface brightness 
elongation to the SW of A1644-N extends towards the main cluster. 




Fig. 3. — Closeup of A1644-S (left) and A1644-N (right) overlayed with isophotal contours of X-ray surface brightness. The physical 
image scale is 0.96 kpc/arcsec. Color gradient is increasing surface brightness from light to dark. 
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some low level of brightness, below which the smoothing 
scale is constant). This method has the effect of apply- 
ing a smaller smoothing scale in areas of bright cluster 
emission and a larger smoothing scale in the low sur- 
face brightness regions, resulting in approximately the 
same relative statistical accuracy across the interesting 
bright regions of the cluster. The surface brightness im- 
ages are then used to fit a mekal model spectrum for each 
image pixel, binned to those same energy bands, and a 
chi square minimization is performed allowing only the 
model temperature and normalization to vary. 

The temperature maps (Figure 7) for the central re- 
gions of the N and S clusters show the cool gas surround- 
ing each cluster's X-ray peak. Their cool cores exhibit 
sharp temperature gradients, with projected gas temper- 
atures going from ~ 1.5 keV near the centers to <~ 6 keV 
outside the surface brightness edge over a radial distance 
of ~ 20" (1" ~ 1 kpc). A small finger-shaped stream of 
warm, less dense (higher specific entropy) gas lies just to 
the west of the cool central peak in A1644-S, extending 
inwards to within ~ 22 kpc of the core, labelled SROI1 
in Figure 4 (left) and Figure 7 (left). In addition to 
cold fronts being created by core gas sloshing, simula- 
tions predict the hotter, outer cluster gas can be drawn 
in towards the center along such finger shaped streams 
(AM06; Poole et al. 2006). We discuss the implications 
of this further in § 5. 

In A1644-N (Figure 7 right), we see a more circu- 
larly symmetric temperature structure near the core, in 
approximate agreement with X-ray brightness contours. 
There is a compression of the X-ray isophotes to the 
southwest and the minimum in the temperature map 
(corresponding to the lowest entropy gas) lies at the in- 
ner edge of this compression. This compression indicates 
that the core in A1644-N is moving to the west with re- 
spect to its surrounding gas. We extract spectra in a 
radial sector across this edge and discuss the results in 
§5.2. 

4.3. Fitting the Surface Brightness Edges 

Figure 5a shows the surface brightness profiles of the 
two sectors in A1644-S (regions SQ2 and SQ3 in Fig- 
ure 4 left ). We observe the edge in each sector, where 
the surface brightness drops abruptly at r^ 12.1" (<~ 12 
kpc) in SQ2 and r~ 32.6" (~ 31 kpc) in SQ3 from their 
respective centers of curvature. We select radial bins in- 
side and outside these edges for spectral extraction, with 
the bin width chosen to maximize signal to noise in each 
region, while maintaining sufficient resolution to resolve 
any radial trends. The fitted gas temperatures (§3.4) are 
shown in Table 1 and plotted in Figure 5b. They con- 
firm that the surface brightness edge in each quadrant 
is indeed a cold front with the deprojected gas tempera- 
ture increasing from kT=1.86^o'2i keV inside the edge to 
kT= 5.74±J || keV in SQ2 and from kT= 1.92±°;^ keV 

to kT= 5.77+^ keV in SQ3. To quantify the magnitude 
of these jumps precisely, we define the temperature jump, 
Jt, as the ratio of the deprojected temperature measured 
from one side of the jump to the other: Jt = T out /Ti n . 
Given that the temperature on either side of the front is 
not uniform (in SQ3 and NQl), we correct for the tem- 
perature gradient by modelling it with a power law on 
each side and obtain its value just inside and just outside 



TABLE 3 
Cold Front Parameters 



Region ID 


Jt 




SQ2 
SQ3 
NQl 


3 nq+ 113 

3 m+ ' 74 
cS - U1 -0.47 

1 SS + ' 66 


2 62+ 1 ' 47 

3 7Q+ 1 ' 32 
1 QS+ 089 



Note. — Table 3: Jumps in tem- 
perature and density across the re- 
gions in Table 1. Ratio of temper- 
ature (col 2) and density (col 3) in- 
side and outside the surface bright- 
ness edges as described in §4.3. Er- 
rors on the density jump are the 90% 
confidence interval. 

the edge radius. We present these in Table 3. 

Since the plasma X-ray emissivity is proportional to 
the square of the gas density (and only weakly depen- 
dent on the temperature for our energy band), the sur- 
face brightness distribution provides a direct estimate 
of the gas density. A discontinuity in the slope of the 
surface brightness profile at the edge therefore implies a 
discontinuity, or a jump, in the gas density. We model 
the density within our spherical sectors encompassing the 
edge with a power law function immediately interior to 
the edge (as in Markevitch et al. 2000 and many later 
works). Outside the edge in A1644-S and A1644-N, a 
second power law profile gives a good fit, whereas for 
the larger deprojection regions a f3 model is used. The 
density models are: 

n e (r)=n 1 J ne ( — ) 71 , 

r 3 

for r < rj in all regions 

for r > r-jin SQ2, SQ3 and NQl 
a 

for r > Tj in S D e P ro 3 and N fle))roj 

where rj is the radius of the edge 4 , 71 is the power law 
slope interior to each edge, 72 is the power law slope 
outside the edge, a is the scale radius of the j3 model, ni, 
712 and 713 combine the respective model normalizations 
with the emission measure normalization for each model 
in their respective regions, and J„ e is the magnitude of 
the density jump. When projected, this density model 
produces an X-ray brightness edge with the characteristic 
profile observed at cold fronts. Inside the edge we limit 
the fit to the immediate vicinity of the edge, as the profile 
near the brightness peak is like to have a more complex 
shape. 

To fit the density change across the jump, this radial 
density profile is squared and integrated along the line 
of sight, multiplied by a factor close to 1 to account for 
the weak dependence of the X-ray emissivity on temper- 
ature, and then fit to the surface brightness distribution. 

4 In general, the centers of the inner and outer models may differ, 
however they are assumed identical 
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Fig. 4. — A1644-S (left) and A1644-N (right) with the spectral extraction regions overlayed. The regions are color and pattern coded to 
match the points plotted in Figures 5 and 6. 
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Fig. 5. — A1644-S: a) Radial surface brightness, b) deprojected temperature, c) electron density, and d) pressure profiles for the regions 
SQ2 (solid x's, medium (red) shading), SQ3 (open dashed squares, dark (blue) shading) in Figure 4 (left), and S£) eproJ (solid squares, light 
(green) shading) in Figure 2 (left). Vertical dashed lines are the radii of the edge in each sector. Spectral extraction and fitting is described 
in § 3.5. Shading in c) and d) are from the 90% confidence intervals on the model parameters. 
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The resulting density fits are shown in Figure 5c and Fig- 
ure 6c along with the model errors. We then compute the 
internal pressure of the gas simply as P = n e T which is 
plotted in Figure 5d and Figure 6d. Figure 5 and Ta- 
ble 3 show that the edge in A1644-S is a very prominent 
cold front, while that in A1644-N is marginally consistent 
with a smooth gas profile. 

5. DISCUSSION 
5.1. Comparison With Simulations 

To assist in interpreting the gas motion which pro- 
duced the spiral morphology we reproduce a snapshot 
from the simulations of AM06 in Figure 9 (also see Fig- 
ure 7 in AM06). Although this simulation is for an R=5 
mass ratio with a gasless subcluster, other simulations 
in AM06 show a qualitatively similar structure on the 
small scales we are comparing here (r < 100 kpc). So 
our conclusions may be equally well drawn from cases 
with initial conditions more appropriate to A1644 (e.g., 
the R = 2, subcluster-with-gas case shown in Fig. 14 of 
AM06; note that on large scales, A1644 does show the 
disturbed morphology consistent with a merger with a 
subcluster with gas). Our choice here was made as it 
affords the highest-resolution morphological comparison 
to already published data. 

At the 1.6 Gyr snapshot, the perturber has already 
passed above the main cluster, from the top right to- 
wards the bottom left in the simulation frames. We find 
the 2.1 Gyr snapshot in the AM06 simulations (700 Myr 
after pericentric passage of the subcluster) to be the 
most similar image to the present state of A1644 (see 
Figure 9 where we show the simulations and our image 
side by side). The cool core gas, which had previously 
been in equilibrium with the cluster potential, was dis- 
turbed by the passage of A1644-N and began sloshing 
around the minimum of A1644-S's gravitational poten- 
tial. The simulations predict that relatively long lived (1 
Gyr) cold fronts can be produced by each passage of the 
gas core about the potential minimum, with their centers 
of curvature at approximately the location of the poten- 
tial minimum of the cluster. Given that we see only one 
such cold front to the east of A1644-S, it appears that 
we are viewing an early stage of gas sloshing. 

5.2. Merger History 

The spiral morphology of core gas sloshing when 
viewed at or near edge-on provides a unique method 
to probe the merger history of the system. Simulations 
(AM06) show that following the spiral structure from the 
outside inwards indicates the direction of the infall of the 
perturbing object. In A1644, the direction of the spiral 
inflow is clockwise (angular momentum vector is into the 
page) and so the perturber must also have a clockwise 
orbit. To produce the observed spiral structure, the per- 
turber most likely passed the main cluster on the east- 
ern side heading in a northerly direction. This is the 
same trajectory of A1644-N as indicated by Reiprich et 
al. (2004). 

Adopting this model, we explore its consequences and 
test whether it agrees with other observational evidence. 
We note that the projected separation of the two clusters 
is <~ 12.4' (~ 700 kpc). The absence of obvious shock 
fronts suggests M< 1.0, so if we assume the velocity of 



the subcluster relative to the main cluster is ~ 900 km 
s" 1 (the sound speed of the 5 keV cluster gas), we may 
estimate a timescale for the pair to arrive at their current 
projected separation. 

Travelling at this average speed, the clusters will have 
obtained this separation ~ 760 / cos i Myr after their 
pericentric passage, where i is the angle of inclination of 
the orbit to the plane of the sky. Returning to the simu- 
lations of AM06, we see that the time from the pericen- 
tric crossing of the two cores to when the first cold front 
morphology is similar to Abell 1644 is ~ 700 Myr (the 
1.4 to 2.1 Gyr snapshots), in rough agreement with the 
timescale for the clusters' current separation. We may 
also compare this to the R=2 and R=5 (subcluster with 
gas) in AM06 (their Figure 9 and Figure 14) and find 
that a similar morphology is produced after ~500 and 
~600 Myr, again in good qualitative agreement with our 
calculation. The isophotes of X-ray surface brightness 
in the northern subcluster, along with its projected tem- 
perature structure suggest that its core is currently mov- 
ing west-southwesterly, back towards A1644-S, although 
they may also indicate merger induced sloshing in that 
core as well. 

5.3. Core Heating 

The X-ray data indicate that the passage of a subclus- 
ter has displaced the cool core of the main cluster from 
the gravitational potential minimum. As it oscillates, the 
newly acquired potentital energy will dissipate into heat. 
It is interesting to compare the rate of such heating with 
the rate of radiative cooling in the core. 

To find the energy necessary to displace its gas core 
from the potential minimum, we require the gravitational 
potential of A1644-S. For an estimate of the gravitational 
potential, we use the mass-temperature relation from 
Vikhlinin et al. (2009) and assume an NFW (Navarro, 
J. et al. 1995) form of the mass profile. Even though 
the cluster is disturbed, this is sufficient for our order 
of magnitude estimate. Obviously, our result is roughly 
dependent on shape, however we note that cold fronts in 
the cluster indicate the presence of a potential cusp (as in 
a NFW profile) as opposed to a shallow cored (/3-model) 
distribution (AM06). 

The potential energy difference between the current 
location of the core (the X-ray brightness peak) and 
the potential minimum (the center of the cD) galaxy is 
U ~ 1.3 x 10 56 erg (see Appendix A for detailed calcula- 
tion). Assuming a fraction / of that energy is converted 
to thermal energy in the gas during the time between 
pericentric crossings of the two cores (~ 2.6 Gyr in the 
AM06 simulations), we may express the resulting heating 
rate as P s i os h = / x 1.5 x 10 39 erg s -1 and compare that 
to the observed X-ray luminosity in the core of A1644-S 
(L core = 2.4 x 10 42 erg s _1 ). Given the strong depen- 
dence of our result on the gas core radius, which is seen 
in projection, this is a lower limit to the available en- 
ergy in the sloshing gas. Additionally, the separation 
between the gas core and the potential minimum is also 
certainly a minimum estimate, as we are not only seeing 
it in projection but we are also assuming that the current 
separation is the maximum separation. 

Core gas sloshing may reheat the core also by drawing 
in higher entropy outer cluster gas along streams towards 
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the center (MV07). A look at the entropy map (Figure 8) 
shows that the stream in A1644-S does contain a higher 
entropy plasma than that found in the core, suggesting 
an origin outside the core. Detailed hydrodynamic sim- 
ulations are needed to estimate the magnitude of this 
effect and will be addressed in a future paper (Zuhone et 
al. in prep). 

6. SUMMARY 

We present an analysis of the 70 ks Chandra ACIS- 
I observation of the double cluster Abell 1644. Based 
on the X-ray and available optical data, we find com- 
pelling evidence that the primary cluster has undergone 
an interaction with its northern companion and that the 
interaction has disturbed the gas around each cluster's 
core. An examination of the X-ray surface brightness 
edge in the southern, main cluster finds that the edge 
is one continuous cold front with a spiral morphology 
arising from the core gas sloshing as seen in the hydro 
simulations of AM06 and across which the temperature 
jumps by a factor of ~3. 

We find that the gravitational potential energy trans- 
ferred to the main cluster's gas core by the perturbing 
subcluster is insufficient to offset the X-ray luminosity. 
Still, the spiral structure in A1644-S presents an interest- 
ing example of how merger-induced sloshing may draw 
higher entropy (i.e. low-density and high T) outer gas 
inside the cooling core, possibly supplying energy to par- 
tially offset cooling. 
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APPENDIX 
APPENDIX A 

Here, we estimate the energy required to displace the core gas in A1644-S to its current (projected) position. The 
general expression for the gravitational potential energy of the gas core is: 



U = GM T (< RjmgR- 1 (Al) 

where R = 1.5" (1.44 kpc) is the distance the core gas has been displaced from the potential minimum, assuming it is 
a solid body for this order of magnitude estimate, Mp(<R) is the total cluster mass within the displaced radius, and 
m g is the mass of the displaced core. 

To find Mt(<R), we fit an NFW profile using the concentration parameter Csoo=3 for a 5 keV cluster (Vikhlinin et 
al. 2006). We then normalize the resultant mass profile at rsoo to the M-T relation in Vikhlinin et al. (2009) to obtain 
M T (<R)=6.7 xlO 9 M . 

For m g , the displaced core gas mass, we integrate the gas density profile obtained from the X-ray surface brightness 
distribution in §4.2 from the center out to the front radius Tj in SQ3 (rj ~ 33" or 31 kpc). From the density profile 
in SQ3 we obtain: 
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Fig. 7. — (upper) Projected temperature maps of the A1644-S (left) and A1644-N (right). Contours are of X-ray surface brightness 
identical to Figure 3. Spectral model parameters for SROI1 are shown in Table 2. Color gradient shown by color bar is temperature in keV. 
(lower) Average of the hi and low 90% confidence interval on projected temperatures. See online version of The Astrophysical Journal for 
color version of this figure. 
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image (0.5-2.5 keV). We see here the finger of gas (SROI1 in Figure 7) is of higher entropy than in the core. Such streams may provide 
a mechanism for increasing core entropy and offsetting cooling. See online version of The Astrophysical Journal for color version of this 
figure. 
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Fig. 9. — Projected temperature (left) and X-ray surface brightness (center) images of the R=5 merger of two clusters reproduced from 
AM06 by permission of the AAS. (right) A smoothed 0.5-2.5 keV image of south cluster of Abell 1644, oriented to match the simulation 
frames (90% clockwise rotation then Y-axis mirror reflection). The physical scales on each of these images are approximately 250 kpc on a 
side. 
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f r 3 m 

m g = Airk / —f-n e (r)r 2 dr (A2) 
Jo 1-1' 

where n e (r) is the same as in §4.3. 

This results in a radial dependence for the core gas mass of: m g (r) oc r\ +l1 . The core gas mass displaced within Tj 
is then: m g ~ 3.1 x 1O 8 M . Substituting M T and m g into Al, we obtain U ~ 1.3 x 10 56 erg which, averaged over the 
time between pericentric core passages (2.6 Gyr) is ~ 1.5 x 10 39 erg s _1 . 



